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ABSTRACT 

Gamma-ray bursts (GRBs) are widely proposed as an effective probe to trace the 
Hubble diagram of the Universe in high redshift range. However, the calibration of 
GRBs is not as easy as that of type-la supernovae (SNe la). Most calibrating methods 
at present make use one or some of the empirical luminosity correlations, e.g., Amati 
relation. One of the underlying assumptions of these calibrating methods is that the 
empirical correlation is universal over all redshifts. In this paper, we check to what 
extent this assumption holds. Assuming that SNe la exactly trace the Hubble diagram 
of the Universe, we re-investigate the Amati relation for low redshift (z < 1.4) and 
high redshift {z > 1.4) GRBs, respectively. It is found that the Amati relation of low-z 
GRBs differs from that of high-z GRBs at more than Str confidence level. This result is 
insensitive to cosmological models. We should be cautious when using Amati relation 
to reconstruct the Hubble diagram of the Universe. 
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1 INTRODUCTION 


Gamma-ray bursts (GRBs) are the most luminous explosions in the Universe since the big bang. The is otropic equiv ¬ 


alent energy they released in a few seconds ca n be as large as 10^® ~ 10®® ergs. For recent reviews, see, e.g., Piran 1 19991 ): 


Meszaxo^200^ 20061 ) : Kumar fc Zhan j ( 20151) . Thanks to their extreme brightness, GRBs are detectable up to redshift z > 9 
(Salv^em 2015). For example, the most distant GRB known today is GRB 090429B, whose redshift is as high as z « 9.4 


i Cucchiara et al.ll201ll ). Due to their high redshift properties, GRBs are often proposed as potential candles to trace the Hubble 


diagram of the Universe in the high redshift range. In fact, GRBs have already been widely used, tog 
such as tvpe-Ia supernovae (SNe la), to constrain the cosmological parameters (ISchaeferl 2003: 

ether with other candles. 

Bloom. Frail & Kulkarni 

2003; 

Xu. Dai & Liand 200,4 Firmani et al.l 200,4 Liang & Zhane 200,4 Firmani et al. 20064 Schaefei 

2007: Liang et al. 

2008; 

Liang & Zhandl2008; IWei & Zhang 2004 IWeilboiC; IWang. Oi & Dai I2OI1I: iGaoozziello et al. 

2012 

: IWei. Wu & Melia 

2013; 

Velten. Montiel & Garneirol 201.4 iGai et al. 2013:lBret6n & Montie]ll2ni.4 Chang: et al. 12014: iGano & .lakobsson l20l4 


luminosities of SNe la make them the ideal distance indicators in tracing the Hubble diagram of the local (low-redshift) 
universe. However, since we have little knowledge about the explosion mechanism of GRBs, the GRB candle is much less 
standard than the SN la candle. 

Nevertheless, one can still calibrate GRBs using the empirical luminosity correlations found in long GRBs. These cor- 


^peak 


relations includes Amati relation (-Ep eak — ^iso) ( Amati et al.l 12003 : Amatil 20031. 200€), Ghir landa relation (-Epeak — E^) 


( Ghirlanda. Ghisellini Lazzati 2004 ). Yonetoku relation {Ep^ak — 7/iso) ( Yonetoku et al. 20041). Li ang-Zhang relation (tf, — 


— Eiso) ( Liang &; Zhang 2005 ). Firmani relation (To.45 — Ep^ak — Liso) ( Firmani et al. 2006bl) . lag-luminosity relation 
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iUss. _ L. iso 1 llNorris. Ma ani fc Boniid] l200o[ ). variability-luminosity relation iV ~ L\so) JPenimore fc Ramirez-Rui3 12OO0I : 


Reichart et al 


200l|), and so on. Among these luminosity correlations, the Amati relation is most widely used. This is partly 


because that the spectrum properties such as the peak energy Speak, the spectrum indices a and /3, and the photon fluence S 
which are necessary to analyze the Amati relation can be easily observed with enough precision, so that the number of available 
GRBs is large. Unfortunately, all of these correlations depend on a specific cosmological model. Therefore, circularity problem 
occurs when using GRBs to constrain the cosmol ogical parameters. R ecently, some model-independent methods were proposed 
to calibrate GRBs, su ch as the Bayesian method ( Firmam_^^_^ 2005), the luminosity distance method ( Ghirlanda et al.ll2004l : 
Liang fc Zhanj 20051 ). and the scatter method ( Ghirlanda et al. 2004 ). etc.. Actually, these methods still can’t completely 


solve the circularity problem. _ 

A completely model-independent method free of circu larity problem is using distance ladder to calibrate GRBs lILiang et al. 
20081 : Liang fc Zhang 20081 : Wei &: Zhanj 20^ : Wei 2O10l ). The main procedures are as follows: Firstly, calculate the distance 


moduli for the low-redshift (e.g., z < 1.4) GRBs by using cubic interpolation from SNe la (e.g., Union2). Thus, the distance as 
well as the isotropic equivalent energy of low -2 GRBs can be derived. Then we can obtain the empirical luminosity correlations 
(such as Amati relation) from the low -2 GRBs. By directly extrapolating the empirical luminosity correlations to high -2 (e.g., 
2 > 1.4) GRBs, we can inversely obtain the distance moduli for high -2 GRBs. Since the distance moduli of SNe la are directly 
extracted from their light curves without involving any cosmologi cal model, this cal ibrating method is of course completely 
model-independent. Rece ntly, a simi lar method was proposed by iLiu fc Weil ll2014h . The only difference is that they used 
the Bade approximation I Fade 18921 ) instead of the cubic interpolation to derive the distance moduli of low -2 GRBs. An 
underlying, but unproven assumption of these methods is that the empirical luminosity correlations are universal over all 
redshifts. If the empirical luminosity correlations evolve with redshift, these calibrating methods might be invalid. In fact, 


Wang. Qi fc Pail 1120111) have already investigated six empirical luminosity correlations in different redshift ranges. They found 


that the slope of Amati relation of high -2 GRBs is smaller than that of low -2 GRBs, although the intercept does not vary 
significantly with redshift. Similar features were found in the rest five luminosity correlations. Due to the large uncertainties, 
they concluded t hat no significant evidence for the redshift evolution of the luminosity correlations was found. However, in 
another paper llLill2007ll . the author has showed that the Amati relation varies with redshifts systematically and significantly. 

In this paper, we focus on checking the redshift dependence of Amati relation. Firstly, we use the Union2.1 dataset 
I Suzuki et al. 20121 ) to constrain the Hubble diagram of the Universe in the redshift range of 2 < 1.4. Then we directly 
extrapolate the Hubble diagram to high redshift range. By assuming that GRBs follow the same Hubble diagram, we can 
calculate their luminosity distance and isotropic equivalent energy. Finall y, we study the Am ati relation for the low -2 (2 < 1.4) 
and high -2 (2 > 1.4) GRBs, respectively. The GRB sample is taken from Liu fc ( 2014 ). This sample consists of 59 low -2 
GRBs and 79 high -2 GRBs. The rest of the paper is arranged as follows: In section we introduce the data and methodology 
that are necessary to our studies. We present our results in section [S] Finally, discussions and conclusions are given in section 

SI 


2 DATA AND METHODOLOGY 


2.1 Union2.1 and Hubble diagram 


SNe la are usually regarded as ideal distance indicators to trace the Hubble diagram of the Universe due to their consistent 
luminosity. The GRB candle, however, is much less standard than the SN la candle, since the explosion mechanism of GRBs 
is still not clearly known. Therefore, SNe la are often us ed to calibrate the distance moduli of GRBs. In this paper, we first 
use the recently published Union2.1 I Suzuki et al. 20121 ) dataset to constrain the cosmological parameters. The Union 2.1 
dataset is a compilation of 580 well-observed SNe la in the redshift range 2 £ [0.015, 1.415]. All the SNe la have high-quality 
light curves, so their distance moduli can be extracted with high precision. The Hubble diagram can be tightly constrained 
by the Union2.1 dataset. 

In the spatially-flat isotropic spacetime, the luminosity distance can be expressed as a function of redshift as 

where c = 3 x 10® m s“^ is the speed of light, Hq — 70 km s“^ Mpc~^ is the Hubble constant at present time, and 

E[z) = H{z)/Ho is the normalized Hubble parameter. In the ACDM model, we have 

E ( 2 ) = 11m(1 + 2 :)® -t- (1 — Hm), (2) 

where Qm is matter density today. In the wGDM model, we have 

E\z) = L!m(1 + zf + (1 - L!m)(1 + (3) 

where w = p/p denotes the equation o f state of dark energy. In the Ghevallier-Polarski-Linder (GPL) parametrization 

i Chevallier fc Polarskil I2OOII : lLinderll2003l ). the equation of state of dark energy is given as Wde = wo + wiz/(I -\- z), and 
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in this case, E{z) can be expressed as i Nesseris fc Perivolaropoulosll2003 : lLazkoz. Nesseris fc Perivolaropoulos 2005) 

E\z) = + zf + (1 - nM)(l + exp • (4) 


In practice, it is convenient to define a dimensionless quantity called distance modulus, that is, 

where “log” is the logarithm of base 10. The best-fit cosmological parameters can be derived by minimize i-®-, 

(Pth “ UohaY 


-E- 

i=l 


j2 


(5) 


( 6 ) 


where pth is the theoretical distance modulus calculated from Eq. fioha is the observed distance modulus, and is the 
measurement error. N = 580 is the number of SNe la in the Union2.1 dataset. 


2.2 GRBs and Amati relation 


The A mati relation is a correlation between i sotropic equivalent energy Ej^o and spec t rum p eak energy in the comoving 
frame Ep^i I Amati et al. 20021 : Amati 20031 . 200dl . It was first discovered by Amati_etah _([2002)hi 12 BeppoSAX GRBs of 
low redshift {z < 2, except one GRB), and was confirmed later in larger samples ( Amatill2003l . 20^ 1. The Amati relation can 
be parameterized as 


, Eiso T , Ep i 

log- = a + 61og 

erg 300 keV 


(7) 


where 

Eiso = 47rd|Sboio(l + z)~^ (8) 

is the isotropic equivalent energy in the 1 keV — 10 MeV energy band, and Sboio is the bolometric fiuence. The uncertainty 
of Eiso propagates from that of Sboio, i.e, 

^Eiso = 47rdicrSboio(l + z)~\ (9) 

The uncertainty from db is absorbed into the intrinsic scatter Uint. Defining 


y = log 


Eiso 

erg 


X = log 


Ep 


300 keV’ 


we can rewrite Eq. 0 as 


The uncertainties of y and x are given as 


y = a + bx. 


1 (^Ei, 


IXE„ 


( 10 ) 


( 11 ) 


( 12 ) 


^ In 10 Eiso ’ In 10 Ep,i ’ 

where “In” represents the natural logarithm. 

The slope and intercept of Amati relation, i.e., b and a in Eq. (HU, can be obtained by directly fitting Eq. HU to the 
observed GRB data. However, the plot of the Amati relation in the {x, y) plane shows significant error bars in both the 
horizontal and vertical axes. Besides, the intrinsic scatter dominates over the measurement errors. Therefore, the ordinary 
least-x^ method does not work. We may get different best-fit parameters depending on whether we minimize the sum of 
squared residuals in the y axis or that in the x axis. To avoid such a problem, we use the fitting method presented in 
D’Agostinil (120051 1. The joint likelihood function for the slope b, intercept a, and intrinsic scatter Uint is given as 


C{aiot,a,b) oc ■ 


X exp 


{yi - a - bxiA 


2(cr: 


-e cr^. I 


A^'Lx. + 

The minus-log-likelihood is given as 

- lnZ:(crint,a, 6) = ^ ^ln(cr?nt + al ) + ^ ^ 2 + const. 


(13) 


(14) 


The sum runs over all the GRB data points. The best-fit parameters are the ones which can minimize the right-hand-side 
of Ea. (ll4ll . Note that the Amati relation depends on the luminosity distance db, which further depends on the specific 
cosmological model. 


The GRB dataset used to our analysis is directly taken fro m (201^. The dataset is a simple summation of 

samples presented in two previously published papers ( Wei 20ld : Qin fc Ghen 2013 1. The two papers collected the data from 
many other published papers. This dataset is at present the largest and most up-to-date sample available to analyze the 
Amati relation. This sample is a collection of GRBs observed by various instruments, such as BATSE, BeppoSAX, HETE-2, 
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Table 1. The best-fit cosmological parameters and their Itr uncertainties from the Union2.1 dataset in three different cosmological 
models. 



Um 

w 

WQ 

Wl 

ACDM 

irCDM 

CPL 

0.2798 ±0.0130 
0.2755 ± 0.0640 
0.2962 ± 0.2224 

-0.9903 ±0.1431 

-1.0090 ±0.2249 

-0.2455 ±2.9514 


Table 2. The intrinsic scatters, intercepts and slopes of Amati relation for low-z and high-z GRBs in three different cosmological models. 




aint 

a 

b 

ACDM: 

low- 2 : 
high-2 

0.3810 ± 0.0433 
0.3133 ±0.0313 

52.7326 ± 0.0568 
52.9459 ± 0.0476 

1.6020 ±0.1004 
1.3008 ±0.1122 

«)CDM: 

low- 2 : 
high-2: 

0.3810 ± 0.0433 
0.3132 ± 0.0313 

52.7330 ±0.0568 
52.9479 ± 0.0476 

1.6023 ±0.1004 
1.3010 ±0.1121 

CPL: 

low- 2 : 
high-2: 

0.3810 ± 0.0433 
0.3135 ±0.0313 

52.7328 ± 0.0568 
52.9423 ± 0.0476 

1.6020 ±0.1004 
1.3005 ±0.1122 


Konus-Wind, Swift, Fermi, and so on. The spectra properties (photon indices, peak energy, fluence, etc.) of these GRBs are 
measnred with enough precision. The redshifts are well determined through the observation of afterglows. The bolometric 
fluence is calculated in the rest frame 1 — 10, 000 keV energy band by using the Band function (Schaefer 2007). Our sample in 
total includes 59 low-z (z < 1.4) and 79 high-z (z > 1.4) GRBs. We analysis the low-z and high-z GRBs separately. To check 
the possible model-dependence of the Amati relation, we investigate it in three different cosmological models, i.e., AGDM 
model, tnCDM model and CPL model. 


3 RESULTS 

To test the Amati relation, the Hubble diagram of the Universe should be known a priori. Thus, we first use the Union2.1 
dataset to constrain the Hubble diagram of the Universe in low-redshift range (z < 1.4). The constraints on the parameters 
of three cosmological models are listed in Table [T] The quoted errors are of Icr. In the fitting procedure, we fix the Hubble 
constant to be Ho = 70.0 km s“^ Mpc“^. Note that the parameters of ACDM model and tnCDM model can be tightly 
constrained. However, the constraint on the CPL model is rather loose, especially on the parameter wi. 

We directly extrapolate the Hubble diagram to the whole redshift range. Assuming that GRBs follow the same Hubble 
diagram, we can calculate the luminosity distance of GRBs from Eq. 0, as well as the isotropic equivalent energy from Eq. ®. 
The Amati relation Eq. m is then used to fit the low-z and high-z GRBs, respectively. The best-fit parameters are listed in 
Tabled The quoted errors are of la. The Upeak — Eiao correlation is plotted in Figure [T] in logarithmic coordinates. Low-z 
and high-z GRBs are denoted by black and red dots with la error bars, respectively. The lines stand for the best-fit results. 
Besides, the la, 2a and 3a contours in the (a, b) plane for low-z (black curves) and high-z (red curves) GRBs are plotted in 
Figured The best-fit central values are denoted by dots. 

From Table [2] and Figure [T] we can see that the slope and intercept of low-z GRBs differ from that of high-z GRBs 
significantly. Low-z GRBs have larger slope, but smaller intercept than high-z GRBs. This can be seen more clearly from the 
contour plot in the (o, 6) plane in Figure [21 There is no overlap between la contours, while only a little overlap between 2a 
contours. The best-fit central values of low-z GRBs do not locate in the 3a contour of high-z GRBs. Similarly, the best-fit 
central values of high-z GRBs locate outside of the 3a contour of low-z GRBs. This means that the Amati relation of low-z 
GRBs differs from that of high-z GRBs at more than 3a confidence level. Another interesting feature is that, there seems 
to be a positive correlation between the slope and intercept for low-z GRBs. While for high-z GRBs, the slope is negatively 
correlated to the intercept. We can also see that the Amati relation is insensitive to cosmological models. 

We note that most high-z GRBs are energetic (Eiso > 10®^ ergs) and have large peak energy {Ep^i > 300 keV). On the 
contrary, low-z GRBs are in general less energetic and have smaller Ep^i than high-z GRBs. Especially, there are two dim 
GRBs with isotropic equivalent energy Eiso < 10®°, while most other low-z GRBs have Eiso > 10®^. The best fit line for 
low-z GRBs is largely affected by the two dim GRBs. GRBs with isotropic equivalent energy small than 10®^ ergs tend to 
be detected only at low redshift, so that the Eiso coverage in redshift may be biased. This may be a reason why the Amati 
relations for low-z and high-z GRBs have very different slopes. 
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log (Sp,i/300 keV) log (Sp,i/300 keV) log (ilp,i/300 keV) 


Figure 1. The Amati relation for low- 2 : (black line) and high -2 (red line) GRBs in different cosmological models. From left to right: 
ACDM model, tuCDM model and GPL model. 



a a a 


Figure 2. The Icr, 2a and 3(7 contours in the (a, b) plane for low -2 (black curves) and high -2 (red curves) GRBs in different cosmological 
models. The central values are denoted by dots. From left to right: ACDM model, ruCDM model and GPL model. 


4 DISCUSSIONS AND CONCLUSIONS 


In this paper, we investigated the Amati relation for low -2 and high -2 GRBs separately in three cosmological models, 
i.e., ACDM model, wCDM model and CPL model. SNe la were nsed to constrain the parameters of background cosmos. It 
was fonnd that the Amati relation of low -2 GRBs differs from that of high -2 GRBs at more than 3cr confidence level. We 
noted that the slope of low -2 GRB is larger tha n that of high -2 GRBs, while high -2 GRBs have larger intercept than low -2 
GRBs. Our results are consistent with that of IWang. Qi fc Pail ( 201 ll ). We also found a positive correlation between the 
slope and intercept for low -2 GRBs. While for high -2 GRBs, the correlation is negative. These features do not depend on the 
cosmological models. We noted that high -2 GRBs are often much energetic than low -2 GRBs. The i?iso coverage in redshift 
may be biased since dim GRBs are only detected at low redshift. 

The Amati relation is often used to calibrate the distance moduli of GRBs. These calibrating methods are based on the 
assumption that the Amati relation is universal over all redshifts. The redshift-dependence of Amati relation, as indicated 
in this paper, puts a challenge on this assumption. However, this does not necessarily mean that the Amati relation indeed 
varies with redshift. Another possibility is that the high -2 Hubble diagram deviates from the low -2 Hubble diagram predicted 
by SNe la. But this is less likely, because the constraints on the Hubble diagram from t he combination o f SNe la and GRBs 
(which are calibrated using the Amati relation) is very close to that from SNe la alone IILiu &: Weill2014l ). 

The Amati relation has intrinsic scatter much larger than the measurement errors. The uncertainties of distanc e moduli of 
GRB s calibrated through Amati relation is about one order of magnitude larger than the uncertainties of SNe la lILiu &: Wei 
201411 . When combining GRBs with SNe la to constrain cosmological parameters, the weight of GRBs is about two orders of 
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magnitude smaller than that of SNe la, since the weight of a data point is inversely proportional to the square of its uncertainty 
(see Eq.®). This is one reason why adding GRBs to the Union2.1 dataset (or any other datasets which h ave much higher 
precision than GRBs) has only small effect on constraining the cosmological parameters IILiu fc Weill2014l ). Another reason 
is that the number of GRBs is much smaller than the number of SNe la. In addition, use GRBs alone to constrain the 
cosmolog ical parameter s will l ead to unreasonable results. All of these arouse us to search for other calibrating methods. For 
example, Basak fc Rao ( 2013l l pointed out that the pulse-wise Amati relation is more robust than the time-averaged one. 
The uncertainties should be much reduced before GRBs can be used, in combination with other candles, to trace the Hubble 
diagram of the Universe. 


ACKNOWLEDGEMENTS 

We are grateful to Y. Sang, P. Wang and D. Zhao for useful discussions. This work has been funded by the National 
Natural Science Fund of China under grants Nos. 11375203, 11305181, 11322545 and 11335012, and by the project of Knowledge 
Innovation Program of Chinese Academy of Sciences. 


References 

Amati L., 2003, Chin. J. Astron. Astrophys. Supp., 3, 455 

Amati L., 2006, MNRAS, 372, 233 

Amati L., et ah, 2002, A&A, 390, 81 

Basak R., Rao A. R., MNRAS, 2013, 436, 3082 

Bloom J. S., Frail D. A., Kulkarni S. R., 2003, ApJ, 594, 674 

Breton N., Montiel A., 2013, Phys. Rev. D, 87, 063527 

Cai R. G., Ma Y. Z., Tang B., Tuo Z. L., 2013, Phys. Rev. D, 87, 123522 

Cano Z., Jakobsson P., 2014, arXiv: 1409.3570 

Capozziello S., Consiglio L., De Laurentis M., De Rosa G., Perna C., Vivolo D., 2012, arXiv: 1206.6700 
Chang Z., Li X., Lin H.-N., Wang S., 2014, Mod. Phys. Lett. A, 29, 1450067 
Chevallier M., Polarski D., 2001, Int. J. Mod. Phys. D, 10, 213 
Cucchiara A., et ah, 2011, ApJ, 736, 7 

Cuzinatto R. R., Medeiros L. G., de Morais E. M., 2014, arXiv: 1412.0145 
D’Agostini G., 2005, arXiv:physics/0511182 [physics.data-an] 

Fenimore E. E., Ramirez-Ruiz E., 2000, arXiv: astro-ph/0004176 

Firmani C., Avila-Reese V., Ghisellini G., Ghirlanda G., 2006a, MNRAS, 372, L28 

Firmani C., Ghisellini G., Ghirlanda G., Avila-Reese V., 2005, MNRAS, 360, LI 

Firmani C., Ghisellini G., Avila-Reese V., Ghirlanda G., 2006b, MNRAS, 370, 185 

Ghirlanda G., Ghisellini G., Lazzati D., 2004, ApJ, 616, 331 

Ghirlanda G., Ghisellini G., Lazzati D., Firmani C., 2004, ApJ, 613, L13 

Kumar P., Zhang B., 2015, Phys. Rep. 561, 1 

Lazkoz R., Nesseris S., Perivolaropoulos L., 2005, J. Cosmol. Astropart. Phys., 11, 010 

Li L. X., 2007, MNRAS, 379, L55 

Li Z., Ding X., Zhu Z. H., 2015, arXiv: 1504.03482 

Liang N., Xiao W. K., Liu Y., Zhang S. N., 2008, ApJ, 685, 354 

Liang E. W., Zhang B., 2005, ApJ, 633, 611 

Liang N., Zhang S. N., 2008, AIP Conf. Proc., 1065, 367 

Linder E. V., 2003, Phys. Rev. Lett., 90, 091301 

Liu J., Wei H., 2014, arXiv: 1410.3960 

Meszaros P., 2002, ARA&A, 40, 137 

Meszaros P., 2006, Rep. Prog. Phys., 69, 2259 

Nesseris S., Perivolaropoulos L., 2004, Phys. Rev. D, 70, 043531 

Norris J. P., Marani G. F., Bonnell J. T., 2000, ApJ, 534, 248 

Pade H., 1892, Ann. Sci. Ecole Norm. Sup., 9, 1 

Piran T., 1999, Phys. Rep., 314, 575 

Qin Y. P., Chen Z. F., 2013, MNRAS, 430, 163 

Reichart D. E., Lamb D. Q., Fenimore E. E., Ramirez-Ruiz E., Cline T. L., Hurley K., 2001, ApJ, 552, 57 
Salvaterra R., 2015, arXiv: 1503.03072 
Schaefer B. E., 2003, ApJ, 583, L67 


© 2015 RAS, MNRAS 000, [IH7] 












Are long gamma-ray bursts standard candles? 7 


Schaefer B. E., 2007, ApJ, 660, 16 
Suzuki N., et al., 2012, ApJ, 746, 85 

Velten H., Montiel A., Carneiro S., 2013, MNRAS, 431, 3301 
Wang F. Y., Dai Z. G., Liang E. W., 2015, arXiv:1504.00735 
Wang F. Y., Qi S., Dai Z. G., 2011, MNRAS, 415, 3423 
Wang J. S, Wang F. Y., 2014, MNRAS, 443, 1680 
Wei H., 2010, J. Gosmol. Astropart. Phys., 08, 020 
Wei J. J, Wu X. F, Melia F, 2013, ApJ, 772, 43 
Wei H., Zhang S. N., 2009, Eur. Phys. J. G, 63, 139 
Xu D., Dai Z. G., Liang E. W., 2005, ApJ, 633, 603 

Yonetoku D., Murakami T., Nakamura T., Yamazaki R., Inoue A. K., loka K., 2004, ApJ, 609, 935 


© 2015 RAS, MNRAS 000, [IHT] 


